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Abstract 



We apply the statistical measure of complexity, introduced by Lopez-Ruiz, Mancini and 
Calbet (LMC), to uniform Fermi systems. We investigate the connection between information 
and complexity measures with the strongly correlated behavior of various Fermi systems 
as nuclear matter, electron gas and liquid helium. We examine the possibility that LMC 
complexity can serve as an index quantifying correlations in the specific system and to which 
extent could be related with experimental quantities. Moreover, we concentrate on thermal 
effects on the complexity of ideal Fermi systems. We find that complexity behaves, both at 
low and high values of temperature, in a similar way as the specific heat. 
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on ! 1 Introduction 

o ■ 

^ ■ Information theory, introduced by Shannon in 1948 to provide answers for some fundamental 
^ . aspects of communications [I], celebrates its first massive application in a quantum system, i.e. 
atoms, almost four decades later [2j. In his seminal work, Gadre stated that information theory 
is a hidden treasure yet to be discovered, which is still true. His prediction has been justified over 
the past twenty years by the extended applications of information-theoretic methods [H El lU El E] , 
in various quantum systems. Information-theoretical methods play an important role, not just 
in the clarification of fundamental concepts of quantum mechanics, but also provide a series of 
results concerning the information content of systems, the presence of interactions, correlation 
with experimental measured quantities, extraction of universal relations etc [21 [71 El El [101 EEH EL21 
fT3l [TH [El [TEl [ITJ, [THl [HI EH EI]- Additionally, many complexity measures have been proposed as 
indicators of complex behavior found in different systems scattered in a broad spectrum of fields 

[22i m E2 esi m eh eh m\ mmmmm^mmmmmmmmmmmmm- 

The statistical measure of complexity Clmc introduced by Lopez- Ruiz, Calbet and Mancini 
(LMC) [22] identifies the entropy or information stored in a system and its distance to the equilib- 
rium probability distribution as the two basic ingredients giving the correct asymptotic properties 
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of a well-behaved measure of complexity. So far several complexity measures have been proposed 
as indicators of complex behavior in various systems, mostly coming from physics, computa- 
tional sciences e.t.c. The LMC measure of statistical complexity is an easily calculable measure 
(compared with other definitions e.g. Kolmogorov's one), defined in the form of the product 
Clmc — SD, combining information S and disequilibrium D. It has indeed the features and 
asymptotical properties that one expects intuitively. It vanishes for the two extreme cases of a 
perfect crystal (perfect order) and ideal gas (perfect disorder). The initial definition of Clmc 
has been slightly modified in a suitable way by Catalan et al [23], leading to the form C = e s D 
applicable to systems described by either discrete or continuous probability distributions. In [23] 
it was shown that the results in both, discrete and continuous cases, are consistent: extreme val- 
ues of C are observed for distributions characterized by a peak superimposed onto a uniform sea. 
Moreover, C should be minimal, when the system reaches equipartition and the minimum value of 
C is attained for rectangular (uniform) density probabilities giving the value C = 1. Additionally, 
C is not an upper bounded function and can become infinitely large. 

LMC complexity is referred to in the literature as shape complexity, since it exhibits larger 
values for complicated patterns of probabilities, as seen, in a first step, intuitively and by inspecting 
the plots. The first investigation of Clmc m quantum many-body systems was carried out in 
atoms, for continuous electron distributions [33] and discrete ones [31]. An alternative definition 
of complexity is the SDL measure T a p [26J (Shiner, Davison, Landsberg), defined and calculated 
in an analogous way as the LMC one. It has been applied in atoms as well, starting from [32] . 
Comments on the validity of SDL and LMC measures are given in detail in Section 4 of [33] , where 
it is stated that a welcome property of a definition of complexity might be the following: If one 
complicates the system by varying some of its parameters, and this leads to an increase of the 
adopted measure of complexity, then one could argue that this measure describes the complexity 
of the system properly. 

The modified version C = e D, also known as shape complexity, satisfies some additional 
and desirable features such as positivity, invariance under translations, rescaling transformations 
and replication. Also, another indication of its internal consistency is the fact that the first two 
g-values of the Renyi entropy [H] are the two defining elements of shape complexity, i.e. e s and 
D [50] . The usefulness of the improved version has been shown in many fields [29j [301 ED E21 [33j 
[Ml EH [361 EZl [381 SHI SQl SH S21 HH S5]. Moreover, the specific measure is suitably tailored 
for quantum systems, described by their very nature probabilistically via density distributions in 
position and momentum spaces, which are necessary for and enable a relatively easy calculation 
of S and D, entering the formulas Clmc = SD or C = e s D. 

The motivation of the present work, in the spirit of the above statements, is to extend our 
previous study of uniform Fermi systems [16j, beyond information entropy, in order to include 
the complexity measure proposed by Lopez- Ruiz et al. [22J, using probability distributions in 
momentum space. In uniform systems the density p = N/V is a constant and the interaction of 
the particles is reflected to the momentum distribution which deviates from the theta function 
form of the ideal Fermi-gas model. Our aim is to connect C, a measure based on a probabilistic 
description and the shape of the corresponding momentum distributions to the phenomenological 
parameters introducing the inter-particle correlations and experimental data (e.g. specific heat). 
It is important to examine how the interaction affects the momentum distribution as well as the 
complexity. An attempt is also made to relate the complexity C with statistical quantities such 
as the temperature. 

The complexity C cannot be measured experimentally, but it is possible, as we demonstrate 
here, to calculate it, starting from a reasonable definition (LMC or SDL) and in our case, by 
employing an information-theoretic method, developed in previous work [IB]. Experiment can 
enter our work through the experimentally measured momentum distribution n(k). Momentum 
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distributions n(k) can be assessed experimentally (e.g. in nuclear matter and liquid helium) via 
deep inelastic scattering at a large momentum transfer. The extraction of n(k) from the measured 
scattering intensity is influenced by the limitations imposed by the experimental resolution and 
the final-state interaction. Hence, the most accurate information on n(k) available to date is likely 
to be one obtained through accurate theoretical calculations, aided by experimental input. 

In the present work, the quantum systems under examination are nuclear matter, electron gas 
and liquid 3 He. The inter-particle interactions of these systems differ in general by many orders 
of magnitude in strength and range. The corresponding potentials, scaled under suitable energy 
and length measures for the different systems, i.e. Fermi energy and inverse Fermi momentum, 
still differ by orders of magnitude. The 3 He system is the most strongly interacting one at short 
distances with an almost-hard-core interaction, while electron gas is the most weakly interacting. 
Nuclear matter lies somewhere in between. Helium and nuclear potentials have relatively weak 
attractive tails. The electronic potential is quite different. It has a weak core (compared with 
3 He and nuclear matter) but its rate of decrease for large r is slow. Thus, at large distances, the 
electronic potential is stronger than the other two. 

Furthermore, the density under investigation affects the assessment of the effect of strong versus 
weak interaction. Characteristic is the example of the electron gas, which is distinguished from 
other systems by the long range nature of the Coulomb interaction. As a result, strong coupling 
prevails in the limit of low density for electron gas, whereas helium and nuclear systems become 
more strongly interacting for higher density regions. In all cases the strength of the interaction 
may be gauged by the depletion of the Fermi sea. Quantitatively, this can be assessed as the 
deviation of Zp from unity, where Zp is the discontinuity gap of the momentum distribution n(k) 
at k = k p , in a uniform Fermi system. The problem of discontinuity of the momentum distribution 
is examined for Fermi-liquids by Migdal [5T]. It is shown that the discontinuity in the momentum 
distribution at k = kp is an inherent consequence of an arbitrary interaction between particles in 
an infinite system. 

The paper is organized as follows. In Section I we present the method leading to the expressions 
of momentum distribution, information content and complexity measures in finite Fermi systems. 
Applications of that expression to nuclear matter, electron gas, and liquid 3 He are made in the 
three subsections of Section II. In the same subsections numerical results are also reported and 
discussed. In Section III the study of the influence of thermal effects on the complexity is made. 
Finally, the concluding remarks and a summary of the present work are given in Section IV. 

2 Momentum distribution, information entropy and com- 
plexity 

The one-body density matrix is the key quantity for the description of the momentum distribution, 
both in infinite and finite quantum systems. It is defined as 



The diagonal elements p(ri,ri) of the density matrix yield the local density distribution, which is 
just a constant p in the case of a uniform infinite system. Homogeneity and isotropy of the system 
require that p(r 1 ,r' 1 ) = p(\ r x — r[ |) = p(r). 

The momentum distribution for fermions is defined by 




ri, r 2 , • • • , r N ) ^(r[, r 2 , . . . , r N ) r 2 . . . r^. 





(2) 
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where v is the single-particle level degeneracy 



v 



2 for electron gas and liquid 3 He 
4 for nuclear matter. 



The normalized momentum distribution, 



J n{k)k = 1, 



is given by the relation 



1 ~/, \ 1 / n-{k), k < k F 



n(k) = —n(k) = — { 1 ' ' V 7 (3) 
1 ' V k K ' V k \ n + (k) } k > kp, 1 ; 



where V k = ~7rk F . The Fermi wave number kp is related with the constant density p = Np 
3/(47rrQ) as follows 

a 2 \ 1/3 / n i \ 1/3 
D7T p \ I 97T 1 x 



In the case of an ideal Fermi gas the momentum distribution has the form 

n (k) = ^-6(k F -k). (5) 

The information entropy in momentum space is given by the relation 

S k = - J n{k)\nn{k)k. (6) 

So, for an ideal Fermi gas, using Eq. fl^D, Sk becomes 

S k = S = \nV k = \n(^^y (7) 

For correlated Fermi systems, the information entropy in momentum space, can be found from 
Eq. (jHJ) if we replace n(k) from Eq. (J3j). Sk is written now [16] 

S k = \nV k -y(jj k 2 h„{k) \nn.{k)dk + k 2 n + {k) \nn + (k)dkj . (8) 

The correlated entropy S k has the form 

Sk = So + S COT , (9) 

where S is the uncorrelated entropy given by Eq. (j7|) and S cor is the contribution of the particles 
correlations to the entropy. That contribution can be found from the expression 

1 — POO 



Scor = — 3 ( / r/"/_(.r)]]i»_(.r)(/.r-i / .r n _(./•) In h + (.r)d.r ) . (10) 
where x = k/kp. 

The disequilibrium .D^ (or information energy, defined by Onicescu [3]), in momentum space 
is defined as 

D fc = / n 2 (k)dk. (11) 



and for an ideal Fermi gas, using Eq. (jSJ), becomes 



D k = D = i-. (12) 
Vh 

In the case of correlated Fermi systems, D k is written as 

D k = ±r^L [ f F k 2 h 2 _{k)dk + k 2 h 2 + {k)dk\ . (13) 

Vk Vk \Jo Jk+ J 

The correlated disequilibrium D k has the form 

D k = D D COT , (14) 
where Do is given in Eq. (I12p and .D cor can be found from the expression 

POO 

„2 2 / \ i , / „2„2 




-Dcor = 3 / x n_(x)<ix+ / x n + (x)dx J . (15) 



The LMC statistical measure of complexity Clmc, in momentum space, is defined as 

Clmc = S k D k , 

where S k is the Shannon information entropy, while D k is the disequilibrium of the system under 
investigation (in momentum space). 

The modified version of the complexity, proposed by Catalan et. al. [23], in momentum space, 
is defined as 

C = H k D k , (16) 
where H k represents the information content of the system defined as 

H k = e s \ (17) 

and ensures positivity of information under any circumstances. 
It is easy to show that 

C = C C cor = e s ™D cor , C = e So D = I. (18) 

The physical meaning of Eq.( fl8l) is clear. In the case of an ideal Fermi gas (see Eq. (JSJ)) C is 
minimal with value Co = 1 (see also |23j). Moreover as pointed out in Ref. [23], C is not an 
upper bounded function and can therefore become infinitely large. From the above analysis it is 
clear that complexity C is an accounter of correlations in an infinite Fermi system. So, the next 
step is to try to find the connection between C and the correlation parameters of the systems. 
The correlations invoke diffusion of the momentum distribution and we expect this effect to be 
reflected on the values of C. 

Here, we apply the low order approximation of the momentum distribution in the case of the 
nuclear matter [521 [53j [M] . For liquid 3 He and electron gas we use the most updated calculations 
for the momentum distribution, that is the results of Moroni et al. [55] and P. Gori-Giorgi et al. 
[56] , respectively. Also, we would like to stress out the fact that our main goal is the accurate 
calculation of the correlated part of information and complexity measures, based in reliable data, 
and not the detailed analysis of the momentum distribution itself. 
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2.1 Nuclear matter 

The model we study is based on the Jastrow ansatz for the ground state wave function of nuclear 
matter 

*(r 1 ,r 2 ,... > r JV )= JJ /(r y ) $(r 1; r 2 , rjv ), (19) 

l<i<j<N 

where =| r — |, $ is a Slater determinant (here, of plane waves with appropriate spin-isospin 
factors, filling the Fermi sea) and f(r) is a state-independent two-body correlation function. Thus, 
the correlation function is taken to be of the Jastrow type [57] 

/(r) = l-exp[-/?V], (20) 

where (3 is the correlation parameter. A cluster expansion for the one-body density matrix p(ri, r^) 
has been derived by Gaudin, Gillespie and Ripka [52j [531 [51] for the Jastrow trial function f[T9|) . 
In the low order approximation the momentum distribution is constructed as 



where 



n hOA (k) = 6(k F — k)[l — k dh + Y(k, 8)] + 8 [k dh Y(k, 2) - [Y(k, 4)] 2 ] , (21) 



P -^+ _ a-k- r k + „ _ r\ k -\ , 

c- l Y{k, /x) = -= + / e"» y + sgn(^) / y, (22) 



and 



^(ff ' ^ t-^. (23) 

while sgn(x) = xj \ x |. The dimensionless Jastrow wound parameter fcdir can serve as a rough 
measure of correlations and the rate of convergence of the cluster expansion is defined as 

k Ail = P |[/(r)-l] 2 r, (24) 

The normalization condition for the momentum distribution is 

J n LOA (k)k 2 k=-k F . (25) 

^From Eq. ( 1241) we obtain the following relation between the wound parameter fcdir and the corre- 
lation parameter j3 

1 fk F x3 



k ^^rAi)- (26) 

It is clear that large values of fcdir imply strong correlations and poor convergence of the cluster 
expansion. In the numerical calculations the correlation parameter (3 is in the interval: 1.01 < 
(3 < 2.482. That range corresponds to 0.3 > fcdir > 0.02 and is reasonable in the case of nuclear 
matter [54]. 

The calculated values of S cor , D cor and C for nuclear matter versus wound parameter fcdir are 
displayed in Fig. 1. S COI and C increase with fc d ir, while D cor decreases. We fitted the numerical 
values of the above quantities, with simple functions of k dir and we find respectively the following 
formulae 

S COT = ak^, a = 2.0586, (3 = 0.6365. (27) 
D cor = 1 + ak^, a = -0.9009, (3 = 0.8325. (28) 
C = 1 + a^ ir e 7fcdir , a = 3.1760, (3 = 0.8257, 7 = -1.6176. (29) 
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The values of the parameters a, (3 and 7, for each case, have been selected by a least squares fit 
(LSF) method. 

Another characteristic quantity, used as a measure of the strength of correlations of the uniform 
Fermi systems, is the discontinuity, Z Fl of the momentum distribution at k/kp = 1. It is defined 

as 

Z F = n(l')-n(l + ) 

The behavior of momentum distribution, as a function of k/kp for various values of the wound 
parameter k^ is indicated in Fig. 2(a). The discontinuity Zp is also displayed in each case. For 
ideal Fermi systems Zp = 1, while for interacting ones Zp < 1. In the limit of very strong 
interaction Zp = 0, there is no discontinuity on the momentum distribution of the system. The 
quantity (1 — Z F ) measures the ability of correlations to deplete the Fermi sea by exciting particles 
from states below it (hole states) to states above it (particle states) [54|. 

The dependence of S cor , D COT and C on the quantity (1 — Zp) is shown in Fig. 3. It is seen 
that S cor and C are increasing functions of (1 — Z F ), while D cor is a decreasing one, as a direct 
consequence of the dependence of the above quantities, on the correlation parameter k^ v . That 
dependence can be reproduced very well by simple expressions as in Eqs. (1271) . (1281) and (1291) 
replacing /c dir by (1 — Zp). For example the expression 

C(Z F ) = l + a(l- Zpf^- Z ^ (30) 

with a = 3.6227, (3 = 0.8024 and 7 = —1.9750 reproduces the numerical values of C very well. 

^From the above analysis we can conclude that LMC complexity C can be employed as a 
measure of the strength of correlations in the same way the wound and the discontinuity parameters 
are used. An explanation of the above behavior of C is the following: The effect of nucleon 
correlations is the departure from the step function form of the momentum distribution (ideal 
Fermi gas) to the one with long tail behavior for k > kp. The diffusion of the distribution leads 
to a decrease of the order of the system (the disequilibrium decreases and the information 
entropy increases accordingly). In total, the contribution of Sk in C dominates over the 
contribution of and thus the complexity increases with the correlations (at least in the region 
under consideration). 

2.2 Electron gas 

We consider as electron gas, a system of fermions interacting via a Coulomb potential. The 
electron gas is a model of conduction electrons in a metal, where the periodic positive potential 
due to the ions is replaced by a uniform charge distribution. The density of the uniform electron 
gas (Jellium) is p = 3/(47rr 3 ) and the momentum distribution is n(x,r s ), a function of both 
x = k/kp and r s = r /aB (where 0,3 = h 2 /me 2 is the Bohr radius). In the Fermi-liquid regime, 
the momentum distribution of the unpolarized uniform electron gas n(x,r s ) is constructed with 
the help of the convex Kulik function G(x) [55] . 

Discontinuity Zp(r s ) is unit for r s = and it is a decreasing function of interaction strength 
r s , as r s increases. The discontinuity gap of the momentum distribution n(k) at the Fermi surface 
narrows as the density decreases, a clear indication that the system becomes more strongly coupled. 
That behavior is due to the fact that the screening of the long-range Coulomb interaction between 
the electrons becomes less effective at lower density. Nuclear matter and atomic 3 He exhibit an 
inverse behavior, where the basic interactions are of short range and Zp decreases as the density 
increases. At large r s , electrons form a Wigner crystal with a smooth n(x, r s ). Interaction strength 
r s <C 1 is the weak-correlation limit and r s 3> 1 is the strong-correlation limit, respectively. For 
intermediate values of r s , a non- Fermi liquid regime may exist with Zp = [56J. 
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We examine the dependence of S cor , D COI and C for the electron gas on the correlation param- 
eter r s , (or p = 3/ (47rr^)) and the discontinuity parameter (1 — Zp). This dependence is displayed 
in Fig. 4. It is seen that, as in the case of nuclear matter, S cor depends on those quantities through 
a two parameter expression of the form 

S CO r(r s ) = ar?, a = 0.1312, (3 = 0.8648. (31) 

The disequilibrium D cor takes the form 

D cor = 1 - ar s , a = 0.0463, (32) 

and the complexity C behaves as 

(7= l + arfe rr ', a = 0.0700, (3 = 1.4144, 7 = -0.1525. (33) 

The most distinctive feature, in the case of electron gas, is the occurrence of a maximum of C for 
high values of the correlations parameter r s in contrast to the case of nuclear matter and liquid 
helium (see bellow), where C is a monotonic increasing function of correlations. For high values 
of r s , the competition between D cor and e Scor (see Eq. f[T8"|) ) leads to the dominance of the trend 
of D cor . More precisely the slope of C is given by 

dC ^ ^ d In D cor + dS COT ^ ^ 

d In Dcor 
dr s 

and (which are always negative and positive respectively when r s increases). It is easy to 
show by applying Eqs. (l3~Tj) and fl32|) (or equivalently, but with less accuracy from Eq. (1331) ) that 
C attains a maximum value C max ~ 1.4052 at r s ~ 9.589. 

The above feature is well reflected also on the 1 — Zp dependence of C as exhibited in Fig. 3. 
(the trend of and D cor is also shown). 

2.2.1 Momentum distribution and complexity for the Wigner Crystal 

In the low-density limit, r s — ► 00, the electron gas undergoes Wigner crystallization. The momen- 
tum distribution of the localized electron is of harmonic-oscillator type and has the form 

/ 1 1 \ 3 ^ 2 - fe2 - f°° 
n(k, r s 00) = I -y I e k F \ Aix I k 2 n(k, r s — > oo)k = 1, (35) 



dr s \ dr s dr s 
Thus, according to Eq. the sign of the slope of C depends on the sum of the terms 



where z = u/k F = 0.24ry 2 [56]. It is easy to prove that after some algebra 



3 3tt 1/2 3 
Scor = - + In — h - In z. (36) 



Additionally, D COT can be found also from the expression 



D cm = ^- 3/2 . (37) 



^From Eqs. (135} . (1371) we find that 



L 



s r , /e\ 3 / 2 

-.Of r ^ (i 



C = e bcoi D COI = [-) ~ 1.5845. (3c 



So, in the low-density limit (very strong correlations) the complexity is independent of the corre- 
lation parameter r s . As we show below the above case is similar to that of an ideal Fermi gas at 
high values of temperature. 



2.3 Liquid 3 He 



The interaction potential for liquid 3 He is very strong at small distances and its core repulsion 
is very hard (but not infinite). Consequently there is a Fermi surface discontinuity of roughly 
Zp ~ 0.3. This small value supports the view that 3 He is the most strongly interacting Fermi 
system of the three systems considered here. The momentum distribution has been calculated 
using diffusion Monte Carlo simulations with the use of trial functions optimized via the Euler- 
Monte Carlo method [55] . 

We examine the dependence of S COT , D cor and C on the density p = 3/(47rr 3 ) and the disconti- 
nuity parameter (1 — Zp), and we present our results in Fig. 5. The dependence of S COT on those 
parameters is described through the following simple two parameter formula (as in the cases of 
electron gas and nuclear matter) 

Scor(po) = «Po, Po = 100p, (39) 

with 

a = 2.2736, = 1.4757, 
while the disequilibrium D cor is described by the function 

OL 

D COT = -. — + 5, a = 0.1321, (3 = 1.6288, 7 = 0.1031, 5 = 0.4062. (40) 

The complexity C behaves as 

C = 1 + ap^e 7p0 , a = 0.0166, (3 = -3.1849, 7 = 5.8556. (41) 

The dependence of S cor , D cor and C on the quantity (1 — Z F ) is seen in Fig. 3. In order to be 
able to compare the results of the various systems, in the case of liquid 3 He the values of S COT 
have been divided by 10 and the values of C by 100. The most distinctive feature of the above 
analysis, in the various systems, is the different behavior exhibited by S cor , D COT and C as function 
of (1 — Zp). For the same values of (1 — Zp) both the values and the trend of these quantities are 
different in those systems. 



3 Thermal effects on complexity in electron gas 

At temperature T = 0, the electrons of the electron gas occupy all the lower available states up 
to a highest one, namely the Fermi level. As the temperature increases the electrons of the gas 
tend to become excited and occupy states of energy of order kT higher than the Fermi energy. In 
general the occupation number of the electron gas n(e), is given by the Fermi-Dirac formula 

n(e) = r— 1 , (42) 



exp 



fc^T ( £ " P) 



where e = (p = hk) is the energy of the electrons, ks is the Boltzmann's constant and /x is the 
chemical potential. For T = 0, the chemical potential of a gas coincides with the Fermi energy 
ep, which is by definition the energy of the highest single-particle level occupied at T = 0. i.e. 

^ = l^^ 2/3 ' (43) 
while the Fermi temperature is defined via the relation 

e F = k B T F . (44) 
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We will examine how information and complexity measures considered in the present study 
are affected by temperature, when it starts to increase above zero. The limits of low temperature 
(quantum regime) and high temperature (classical regime) are studied separately in the following 
subsections. 

3.1 Quantum regime (T < Tp) 

With the term low energy we refer to the limit T <C Tp, since there is only a single characteristic 
value of temperature, the Fermi temperature. In a first approximation, the chemical potential for 
that limit, is [581 EH [60] 



fi = ep 



12 VTf 



(45) 



and so Eq. (j4"2j) becomes 



< x ) = r— — 7—1 > ( 46 ) 



exp 



i (^-i + ge) 



+ 1 



where x = (e/ep) 1 / 2 = k/k F , £ = T/T F < 1 and J °° x 2 n(x) x = 1/3. 

Following the same procedure as in Section II, the information entropy Sk of the electron gas 
at temperature T <C Tp is written 

Sk = Sq + Sthermai; (47) 

where S'o is given by Eq.© and 



^thermal = -3 / X 2 n(x) In Tl(x) X. (48) 
J 

In a similar way we find that Dk is written as 

D k = A,Ahermal, (49) 

where D is given in Eq. (fi~2l) and .Dthermai can be calculated also using the expression 

POO 

Ahcrmal = 3 / X 2 ^ 2 (x)x. (50) 



./0 

Finally, it is easy to show that 

C = CoC cor = e thermal ^thermal j C*0 = 1- (51) 

It is worthwhile to notice that the correlations between the Fermi particles invoke a disconti- 
nuity to the momentum distribution at k = k F (see Fig. 2(a)), while the thermal effect causes just 
a slight deviation from the sharp step function form at T = 0. This is shown in Fig. 2(b), where 
the momentum distribution for an ideal electron gas, at various values of T/Tp has been plotted 
versus kjkp. The origin of the two effects (correlations and temperature) is different and it is seen 
that they influence in a different way the momentum distribution. So, it may be appropriate to 
study qualitatively and also quantitatively the above effects on the various information measures 
and complexity. 

The calculated values of Sthermai for various values of the temperature in the low energy limit 
(T ^ Tp) are shown in Fig. 6. It is seen that ^thermal is an increasing linear function of the 
temperature. The linear equation 

T 

Sthermai = a « = 2.5466 (52) 
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reproduces very well the calculated values of S thermal- That expression of the information entropy 
is similar to the expression which gives the thermodynamical entropy, Ste, for T <C Tp. Ste in 
the low temperature limit has the form [581 [60] 



7T 2 T 

Ste = ^- Nk B — . (53) 
2 1 p 

^From Eqs. (|52|) and (|53j) . we can find a relation between the two entropies in the limit T <C Tp. 
The corresponding relation is of the form 

^thermal = T "TTJ — • (54) 

7i z IN Kb 

The calculated values of D thermal, displayed in Fig. 6, are well reproduced by the formula 



Ahermal = e 1M7 \ T J>). (55) 

Finally, the calculated values of C as a function of T/Tp are shown in Fig. 6. C is an increasing 
function of T/Tp and exhibits a linear trend which is fairly reproduced by the relation 

T 

C= 1 + 0.962—. (56) 
Tp 

3.2 Classical regime ( T > T F ) 

In the classical case (where the density is low and/or the temperature is high and it is assumed 
that n(k) <C 1) relations can also be established between the information measures and complexity 
with temperature. In that case the momentum distribution has the gaussian form [59] 



a\3/ 2 „, 2 n 



2 



n(jfe) = ( - ) ' e~ ak , a = — — -, (57) 

\7T/ 2mKBT 

and is normalized as J n(k) k = 1. The above expression is also written as 

/ 1 i \ 3 / 2 _t 2 i 
"(*)= 7-TT e ^ £ = T /Tp. (5* 



^From Eqs. (j4"8~l) and ( |58j) . the following relation, connecting Sthermai with the temperature, can 
be found [T6] 

thermal = | + In — h | In — . (59) 

z 4 Zip 

Additionally, we can also obtain -Dthermai from the expression 

ol/2 / J> \ ~ 3 / 2 

^thermal 

/From Eqs. §9J, (J6HD we find that 



Ahermal - ~^T/2 if ' ^ 



/ 6 \ 3/2 

C = Cthermal = e St " Ahermal = (-) ^ 1-5845. (61) 

The physical meaning of Eq. (ISTl) is very clear. For high values of temperature (T ^> Tp), 
where the momentum distribution of the gas is described by a Gaussian function, the complexity 
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C is independent of T and takes a constant value. Our finding show that the case of an ideal Fermi 
gas at T ^> Tp, is in contrast to the case of correlated Fermi gas at T = 0, where, in general, C 
is not an upper bounded function. Moreover, the complexity has different trend in the quantum 
mechanical limit (T — > 0) compared to the classical limit (T — > oo). In the classical limit, C is not 
affected by the temperature variation and is a constant of the system. However, for low values of 
T, complexity exhibits strong temperature dependence. 

In Fig. 6 we plot the complexity versus £ = T /Tp both for low and high values of temperature. 
Actually, we calculate the chemical potential for each value of £ and consequently, we know exactly 
the momentum dependence of the occupation numbers given in Eq. ( 1421) . The results confirm the 
numerical approximation of C given in Eq. fl56|) for low values of £ as well as the analytical 
prediction of Eq. (I6T1) for high values of £. 

It is worthwhile to notice that the temperature dependence of C is similar to that of the specific 
heat Cy in an ideal Fermi gas [60]. More precisely, Cy is a linear function of T for T <Tp, while 
it approaches 3/2Nk B as T — > oo [60]. In order to illustrate the above result we display in Fig. 7 
the dependence of C and the specific heat Cy on T (in units Tp). This behavior can be explained 
as follows: Let us consider two momentum distributions n(k), one for T = and another one 
for T > 0. They are in essence different, because for T > a certain number of fermions are 
excited above the Fermi level ep. Specifically, fermions with energies of the order of k B T below 
€f are excited to energies of the order of k B T above ep. However, this holds only for fermions 
with energies within about k B T of the Fermi level, while those with other values of energy have 
no place to go -the states are occupied [6T] . 

To sum up, thermal effects lead to a blurring of the Fermi surface that means the distribution 
function n(e) drops over the range ep ± k B T. As a consequence, in the case of low temperature 
limit, only the momentum distribution (or the occupation number) close to the Fermi surface is 
affected by T and this leads to a linear dependence of C on T. For T ^> Tp, the momentum 
distribution is affected for both low and high values of k in such a way so that the complexity C 
tends to a constant. 

4 Conclusions 

In the present work the recently proposed statistical measure of complexity C has been an issue 
under consideration. Information theoretical measures (information entropy and disequilibrium) 
as well as a statistical measure of complexity C have been calculated, in momentum space, for 
realistic Fermi systems, that is nuclear matter, electron gas and liquid helium. The dependence 
of the above measures on the strength of the correlations has been analyzed and displayed. From 
the present analysis it should become clear that the values of the above quantities could be used 
as a measure of the particle correlations of Fermi systems. We have found that the complexity is 
an increasing function of the correlations both for nuclear matter and liquid helium as expected 
intuitively. However, in the case of electron gas, complexity exhibits a different slope and in 
fact, the function of C(r s ) has a maximum for a specific value of the correlation parameter r s . 
Additionally, we have found the interesting result that for very strong correlations, where electron 
gas undergoes Wigner crystallization, the complexity is independent of the correlations and takes 
a constant value. In order to have a common measure of the various information properties, we 
have displayed them as a function of the discontinuity gap, (1 — Zp). The most distinctive feature 
of the above analysis, in the various systems, is the different behavior exhibited by S COT , D cor and 
C as functions of (1 — Zp). For the same values of (1 — Zp) both the values and the trend of these 
quantities are different in the various systems. Considering that, under certain circumstances, Zp 
can be estimated experimentally, we obtain a first indication that information properties may be 
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related with experimental results. 

Temperature also affects the momentum distribution of an ideal Fermi gas and consequently 
all the related information properties. We have found that for low values of the ratio T/Tp the 
complexity is a linear function of T. However, in the high temperature limit (the well known 
classical Maxwell-Boltzmann distribution), complexity is independent of T and takes a fixed value 
(exactly the same as in the case of Wigner crystallization). Thus, regardless of the reason that 
causes the momentum distribution to exhibit gaussian type dependence on the momenta k, the 
value of the complexity is constant. Furthermore, we have seen that the temperature dependence 
of C is similar to that of the specific heat Cy in an ideal Fermi-gas [60J, both for low and high 
values of T. This is a second indication that one can relate the statistical measure of complexity 
C with experimental data (as the specific heat Cy). However, further work is called for before we 
establish a clear connection between information theoretical measures and experimental data. 

As an epilogue, we would like to mention that physicists have been carrying out research for 
decades, going beyond the mean-field description of quantum many-body systems, by taking into 
account correlations among particles, a very important factor indeed towards a better under- 
standing of these systems. The effect of correlations is connected intuitively with the concept of 
complexity, in a qualitative, and somehow vague way. The present work contributes to a quantifi- 
cation of complexity C in correlated Fermi systems, based on previous research for the information 
entropy of the same systems [16]. It turns out that C(T) and the specific heat Cy(T) are similar 
functions of the temperature T as seen in Fig 7(a). In Fig. 7(b) we plot Cy{C). The dependence 
of Cy on C is approximately linear. In fact there appear two regions of linear dependence with a 
different slope, separated by a cross. The fitted expressions are Cy = —1.7353 + 1.7114C (region 
A) and Cy = —6.3777 + 5. 0783C (region B). In a sense, one may state that Cy{T) can serve as an 
index reflecting the expected increase of complexity as T increases. Furthermore, it is seen that 
for temperatures T ^> Tp complexity reaches a plateau (saturation) i.e. it can no longer increase. 
The textbook definition of the specific heat is the measure of heat energy required to increase the 
temperature of a unit of quantity of a substance by a unit degree. It is noted that here, we observe 
an empirical connection of such an " energy-like" quantity with complexity C calculated employing 
information entropy, which is not related directly to the energy of the system in contrast to the 
traditional concept of thermodynamic entropy. 
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Figure 1: S COT , D COI and C of nuclear matter versus the correlation parameter fcdir- The lines 
correspond to the expressions (1271) . (12H1) . (12"9|) . with the parameters derived by the least squares 
fit method. 
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Figure 2: (a) The momentum distribution for correlated nuclear matter versus k/kp for various 
values of the correlation parameter k^ r (b) The momentum distribution of an ideal electron gas 
versus k/kp for various values of the ratio T/Tp. 
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Figure 3: (a) S cor (b) D COT and (c) C for nuclear matter, electron gas and liquid helium versus the 
discontinuity parameter (1 — Zp). In the case of liquid helium, the values of S COI are divided by 
10 and the values of C by 100. 
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Figure 4: S COI , D cor and C of electron gas the correlation parameter r s . The lines correspond to 
the expressions (I3TI) . (|32|) . (1331) . with the parameters derived by the least squares fit method. 
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Figure 5: S cov , D COI and C of liquid helium versus the correlation parameter p . The values of S COI 
are divided by 10 and the values of C by 100. The lines correspond to the expressions ( |39i) . (HP]) . 
( I4T1) . with the parameters derived by the least squares fit method. 
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Figure 6: Sthermai, -Dthermai and C of an ideal electron gas versus the ratio T/Tp. 
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Figure 7: (a) The complexity C and the specific heat Cy (in units fceiV) of an ideal electron gas 
versus the ratio T/Tp (b) The specific heat Cy versus complexity C. 
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